Numerical Continuation of Bifurcations of Limit Cycles in MATLAB

cl matcont and matcont are matlab continuation packages for the interactive numerical study of a range of parameterized nonlinear dynamical systems, in particular ODEs. matcont is an interactive graphical package and cl matcont is a commandline version. Both packages allow to compute curves of equilibria, limit points, Hopf points, limit cycles, flip, fold and torus bifurcation points of limit cycles. We discuss computational details of the continuation of limit cycles and flip, fold and torus bifurcations of limit cycles in matcont and cl matcont using orthogonal collocation. Instead of the more commonly used fully extended systems we use minimally extended systems. We further describe the use of the matlab sparse matrix routines and the initialization and adaptation of the bordering vectors that are essential in minimally extended system. Finally, we compare the use of the minimally and the fully extended systems in the matlab environment.